Integrated processing methodology to identify the spatial variability of accumulated infiltration for agricultural soil

The infiltration of water in the soil, and its variation in space, is essential to establish the irrigation schedule for crops and to evaluate the possible degrading effects on the soil. The objective was to develop an integrated processing methodology in Rstudio to identify the spatial variability of the accumulated infiltration, in two phases related to pea crops. Field sampling was carried out on a rectangular mesh with 48 points per moment, using double infiltrometer rings. The data were evaluated by means of geostatistical tools adjusted with programming code in Rstudio, defining the relationships between the magnitudes of the accumulated infiltration, for different test instants, without the need to make statistical adjustments to the normality of variables, discriminated over a period between 1 and 80 minutes. The results suggest the existence of spatial variability of the accumulated infiltration in the two evaluated phases, considering that most of the analyzed data were adjusted to multiple variance models, maintaining a degree of spatial dependence, and validating the effectiveness of the adjusted methodology developed and implemented. The spatial relationships were corroborated by means of contour maps, where the spatial variation of the accumulated infiltration between the two identified cultivation moments was observed. The reliability of the interpolation by the Ordinary Kriging method was verified by generating variance maps, establishing the degree of homogeneity of the interpolation. The variability of infiltration confirms the validity of the adjusted methodology implemented.

La infiltración del agua en el suelo y su variación espacial es fundamental para establecer la programación de riego en los cultivos y evaluar los posibles efectos degradativos en el suelo. El objetivo fue desarrollar una metodología de procesamiento integrado en Rstudio, para identificar la variabilidad espacial de la infiltración acumulada, en dos fases para un cultivo de arveja. El muestreo de campo se adelantó sobre una malla rectangular georreferenciada con 48 puntos, por cada momento, utilizando anillos infiltrómetros dobles. Los datos fueron evaluados por medio de herramientas geoestadísticas, ajustadas con código de programación en Rstudio, definiendo las relaciones entre las magnitudes de la infiltración acumulada, para diferentes instantes de prueba, sin la necesidad de realizar ajustes estadísticos de normalidad de variables, discriminados en un periodo entre 1 y 80 minutos. Los resultados sugieren la existencia de variabilidad espacial de la infiltración acumulada en las dos fases evaluadas, considerando que, la mayoría de los datos analizados, se ajustaron a múltiples modelos de semivarianza, manteniendo grados de dependencia espacial, particularmente, respecto al máximo valor acumulado de infiltración, validando la eficacia de la metodología ajustada. Las relaciones espaciales fueron corroboradas con mapas de contorno, en donde se observó la variación espacial de la infiltración acumulada entre los momentos de cultivo identificados. La confiabilidad de la interpolación por el método Kriging ordinario, se verificó mediante la generación mapas de varianza, estableciendo el grado de homogeneidad de la interpolación. La variabilidad de la infiltración confirma la validez de la metodología ajustada implementada.

Spatial variability of soil fertility under agroforestry system and native forest in eastern Amazonia, Brazil

The usage of spatial tools might be helpful in the optimization of decision-making regarding soil management, with technologies that assist in the interpretation of information related to soil fertility. Therefore, the present study evaluated the spatial variability of chemical attributes of the soil under an agroforestry system compared to a native forest in the municipality of Tomé-açu, Eastern Amazon, Brazil. Soil samples were performed at 36 points arranged in a 55 x 55 m grid. The soils were prepared and submitted to analysis in order to determine pH in H2O, exchangeable calcium, magnesium, potassium and aluminium, available phosphorus, potential acidity, organic matter, bases saturation and aluminium saturation. For each soil attribute, the spherical, gaussian and exponential models were adjusted. After the semivariograms fitting, data interpolation for assessment of spatial variability of the variables was performed through ordinary kriging. The spherical and gaussian models were the most efficient models in estimation of soil attributes spatial variability, in most cases. Most of variables presented a regular spatial variability in their respective kriging maps, with some exceptions. In general, the kriging maps can be used, and we can take them as logistical maps for management and intervention practices in order to improve the soil fertility in the study areas. The results principal components indicate the need for integrated management of soil chemical attributes, with localized application of acidity correctors, fertilizers and other types of incomes, using the spatial variability of these fertility variables.

Comportamiento espacial de indicadores fisicoquímicos y calidad del suelo en plantación de Theobroma cacao en Padre Abad, Ucayali, Perú

El cultivo de cacao en el trópico peruano presenta bajos rendimientos, lo que exige buscar alternativas para aumentar la producción. Frente a esta situación, los estudios geoespaciales facilitan el diagnóstico de fertilidad y la aplicación de planes de fertilización. Por lo tanto, se evaluó la distribución espacial de indicadores fisicoquímicos y calidad del suelo en Padre Abad, región Ucayali, Perú. Se evaluaron indicadores fisicoquímicos y la calidad del suelo y se sometieron al análisis geoestadístico, a través del interpolador Kriging, encontrándose coeficientes de variación bajos para arcilla, pH y potasio K y medios para arena, limo, materia orgánica (MO), nitrógeno (N), fósforo (P), potasio (K+), calcio (Ca2+), magnesio (Mg2+), aluminio (Al3+), capacidad de intercambio de cationes (CIC), saturación de bases (SB), acidez cambiable (AC) y saturación de aluminio (SAl). La calidad del suelo varió entre baja a muy baja, con una distribución espacial de 52,24 y 47,76 %. El análisis de componentes principales encontró significancia para arena, limo, MO, N, K+, Al3+, CIC, %AC, %BC y %SAl, todos presentan variación espacial, según su nivel de fertilidad, excepto K+. Los modelos de interpolación con mejores ajustes fue el exponencial y lineal con dependencia espacial fuerte para arena, limo y K, moderada para MO, N, Al, CIC, %SAl e índice de calidad del suelo y débil para %BC y %AC, con una eficiencia de predicción confiable. Se encontró alta variabilidad espacial con valores medios de los diferentes indicadores de fertilidad bajos, no favorables para el desarrollo del cultivo de cacao.

Cocoa cultivation in the Peruvian tropics has low yields, which requires seeking alternatives to increase the production. Faced with this situation, geospatial studies facilitate the diagnosis of fertility and the efficient application of fertilization plans. Therefore, the spatial distribution of physicochemical indicators and soil quality in Padre Abad, Ucayali region, Peru were evaluated. Physicochemical indicators and soil quality were evaluated and subjected to geostatistical analysis through the Kriging interpolator, finding low coefficients of variation for clay, pH and potassium (K+), and medium for sand, silt, organic matter (OM), nitrogen (N), phosphorus (P), potassium (K+), calcium (Ca2+), magnesium (Mg2+), aluminum (Al3+), cation exchange capacity (CEC), base saturation (BS), exchangeable acid (EA) and aluminum saturation (AlS). Soil quality ranged from low to very low quality with a spatial distribution of 52.24 % and 47.76 % respectively. Principal component analysis found significance for sand, silt, OM, N, K+, Al3+, CEC, %EC, %BS and %AlS, all showing spatial variation according to their fertility level, except K+. The interpolation models with best fits were the exponential and linear with strong spatial dependence for sand, silt and K, moderate for MO, N, Al, CEC, %AlS and soil quality index (SQI), and weak for %BS and %EC, with reliable prediction efficiency. The research found high spatial variability with low mean values of the different fertility indicators, not favorable for cocoa crop development.

Influence of spatial variability of soil chemical attributes on the nutritional status and growth of the rubber tree

The productive potential of the rubber tree (Hevea brasiliensis) is dependent on its genetic composition, in addition to edaphoclimatic factors and management practices. However, as soil properties are not homogenous, knowing the spatial variability of soil attributes would be important to increase productivity and reduce production costs. In this context, the objective of this study was to determine the spatial variability of chemical attributes of the soil and its influence on the nutritional status and growth of rubber tree clones. Clones FX 3864, FDR 5788, CDC 312, and RRIM 600 were planted at Jaturnaíba Farm, in the municipality of Silva Jardim, Rio de Janeiro State, Brazil. The sampling sites were distributed at a spacing of 20 × 20 m on the northern and southern sides of the relief. The chemical attributes of the soil (pH, Ca2+, Mg2+, K+, P, Al3+, H+Al, sum of bases, cation exchange capacity, and base saturation) were evaluated at a depth of 0­20 cm in the different clone plantations. Additionally, the N, P, K, Ca, and Mg content as well as trunk circumference and total plant height, were also evaluated. Geostatistics was used to determine the spatial variability of the soil and clone attributes, while Ordinary Kriging was used to draw variability maps of the variables. A difference in the distribution of the variables, which was dependent on the slope of the relief, was detected through the maps. The southern side presented better conditions as some degradation was observed on the northern side. Certain soil characteristics influenced the distribution of the attributes of the planted clones; for example, the low concentration of Ca2+ in the soil caused Ca deficiency in the FX clone on the southern slope, indicating that liming did not supply enough nutrients for this clone. Our results showed that the variability in soil attributes influenced the nutritional status and growth of the rubber tree clones, indicating that variability maps can guide the planting and management of the rubber tree, providing more efficient management.

Carnauba (Copernicia prunifera) straw as an alternative bedding material for dairy cows housed in a Compost barn system / [Bagana de carnaúba como material de cama alternativo para vacas leiteiras alojadas em sistema Compost barn]

The replacement of bedding in compost dairy barns (CB) comprises a recurrent management practice, but bedding materials are often not readily available in all regions and the choice of alternative materials is necessary. The objective was to evaluate the thermal attributes of carnauba straw (CS) bedding in compost dairy barn facilities. Environmental monitoring operations were performed at a commercial farm located in Northeast Brazil. Mini weather stations were used to evaluate environmental variables. The THI was evaluated as one of the comfort parameters. Analysis of the spatial distribution of bed surface temperature (BST) in the CB was performed using geostatistical techniques. The cows remained out of the comfort zone according to THI results. The BST indicated satisfactory performance and from the thermal point of view can be used as alternative bedding material in CB facilities. However, it was observed that the CS showed fast biomass degradation compared to conventional materials, widely known. In addition, inadequate temperature values (< 45°C) were found in the deeper of the CS bed, signaling higher risks of pathogenic microbial activity. Additional studies are needed for searching the proper management plans that increase the life span of the bed formed by carnauba straw.(AU)

A reposição de cama na instalação Compost Barn (CB) compreende uma prática recorrente de manejo, porém os materiais de cama muitas vezes não estão facilmente disponíveis em todas as regiões, sendo necessária a adoção de materiais alternativos. Objetivou-se avaliar os atributos térmicos da cama à base de bagana de carnaúba (BC) para instalações CB. Foram realizadas operações de monitoramento ambiental em uma fazenda comercial situada no nordeste do Brasil. Miniestações meteorológicas foram utilizadas para avaliação de variáveis ambientais. O ITU foi avaliado como um dos parâmetros de conforto. A distribuição espacial da temperatura superficial da cama (TSC) no galpão foi analisada usando técnicas da geoestatística. As vacas mantiveram-se fora da zona de conforto, conforme resultados do ITU. A TSC indicou desempenho satisfatório e, do ponto de vista térmico, pode ser utilizada como material alternativo para cama em instalações CB. Entretanto, observou-se que a BC apresentou rápida degradação da biomassa em comparação a materiais convencionais, amplamente conhecidos. Além disso, valores inadequados de temperatura (< 45°C) foram encontrados nas camadas profundas da cama de BC, sinalizando maiores riscos de atividade microbiana patogênica. Estudos adicionais tornam-se necessários para a busca de planos de manejo que aumentem o tempo de vida útil da cama formada por bagana de carnaúba.(AU)

Selection of Interpolators to Predict Populations of Tibraca limbativentris in Irrigated Rice

Abstract Kriging is a method that estimates values in places not sampled from different interpolators, therefore, widely used to predict the spatial distribution of organisms. However, the different interpolators may vary in performance depending on the organism under study or the area evaluated. The aimed study to compare the ordinary kriging and inverse of distance weighted interpolation methods, applied to the spatial distribution of population density of Tibraca limbativentris in irrigated rice. This study was carried out in Santa Maria, RS, Brazil, in two fields with areas of 1.3 ha and 6.2 ha, respectively. Seven evaluations of the population density of T. limbativentris were carried out, corresponding to the period from sowing to maturation. In these areas the adults of T. limbativentris were quantified and the sum used for the statistical and geostatistical analysis. The sample population of T. limbativentris was submitted to different semivariograms, which were selected through cross-validation. The sample population of T. limbativentris was submitted to different semivariograms, selected by means of cross-validation. Once selected, semivariograms were used in both tested interpolation methods. From the results it was concluded that the ordinary kriging interpolation method performed better in all evaluations performed in both areas. Therefore, we recommend its use for estimating the population density and spatial distribution of T. limbativentris in the irrigated rice throughout the crop phenology. Using appropriate interpolation methods, localized management can be used, reducing costs for controlling this pest and increasing the sustainability of the environment.

Efficiency of cooling systems in broiler houses during hot days / Eficiência dos sistemas de resfriamento em aviários durante dias quentes

ABSTRACT: This study evaluated: (1) the vulnerability of broiler houses with different cooling systems, and (2) the spatial distribution of environmental variables during hot days. Four potentially vulnerable commercial broiler houses in southern Brazil were selected according to the following parameters: absence or presence of different cooling systems, broilers older than 28 days, and outside air dry-bulb temperature over 30°C. Broiler house vulnerability was classified according to the cooling and mechanical ventilation system: cellulose pad cooling (CPC), sprinkling (SPK), fogging (FOG), and mechanical ventilation without evaporative cooling system (VTL). The air dry-bulb temperature (Tdb, °C) and relative humidity (RH, %) were recorded every 10 min. For each broiler house, we evaluated: (1) relative cooling efficiency (RCE) and (2) inside spatial distribution of microclimate variables using a geostatistical technique. The CPC and SPK did not differ (P<0.05) in RCE (81.6% and 80.7%, respectively), but both differed from FOG (23.8%) and VLT (1.87%) systems. The highest variations in indoor Tdb were recorded in the FOG (7 °C), followed by the SPK (4 °C) and CPC (3 °C). In the CPC, there was an increase in RH from the middle to the end of the broiler house near the exhaust fans. In conclusion, the relative cooling efficiency and the inside spatial distributions of environmental variables in the broiler houses were influenced by the existing cooling system.

RESUMO: Os objetivos deste estudo foram (1) comparar a vulnerabilidade de aviários com diferentes sistemas de resfriamento, e (2) avaliar a distribuição espacial interna das variáveis ambientais durante dias quentes. Quatro aviários comerciais potencialmente vulneráveis no sul do Brasil foram selecionados de acordo com os seguintes parâmetros de elegibilidade: ausência ou presença de diferentes sistemas de resfriamento; frangos de corte com mais de 28 dias de idade; e temperatura de bulbo seco do ar externa acima de 30 °C. A vulnerabilidade do aviário foi classificada de acordo com o sistema de resfriamento e ventilação mecânica, sendo: resfriamento por pad cooling (CPC), aspersão (SPK), nebulização (FOG) e ventilação mecânica sem sistema de resfriamento evaporativo (VTL). A temperatura de bulbo seco do ar (Tbs, °C) e a umidade relativa do ar (UR, %) foram coletadas por registradores de dados autônomos a cada 10 min. Para cada aviário foram avaliados: (1) eficiência relativa de resfriamento (RCE) e (2) distribuição espacial interna das variáveis ambientais por técnica de geoestatística. O CPC e SPK não diferiram (P>0,05) na RCE (81,6% e 80,7% respectivamente), mas ambos diferiram do FOG (23,8%) e VLT (1,87%). As maiores variações na Tbs do ar interno foram registradas no FOG (7 °C), seguido pelo SPK (4 °C) e CPC (3 °C). No CPC, houve um aumento da umidade relativa do meio para o final do aviário, próximo aos exaustores. Em conclusão, a eficiência relativa de resfriamento, assim como a distribuição espacial das variáveis ambientais internas foram influenciadas pelo sistema de resfriamento de cada aviário.

Distribución espacial de antracnosis (Colletotrichum gloeosporioides Penz) en aguacate en el Estado de México, México / Spatial distribution of anthracnose (Colletotrichum gloeosporioides Penz) in avocado in the State of Mexico, Mexico

Resumen El aguacate (Persea americana) es una especie cuyo cultivo es de gran importancia nutricional y económica para México; sin embargo, como cualquier otro cultivo, a menudo se ve afectado por plagas y enfermedades que limitan su comercialización a nivel mundial. El hongo fitopatógeno Colletotrichum gloeosporioides es el agente causal de la antracnosis en el aguacate y se manifiesta en las etapas tempranas del desarrollo del fruto, así como en poscosecha y durante el almacenamiento, en condiciones de alta humedad relativa (80%) y temperaturas desde los 20 ◦C. Las pérdidas económicas a causa de este hongo pueden ser de hasta el 20% de la producción. En el presente estudio se aplicaron métodos geoestadísticos para definir la distribución espacial de antracnosis en frutos de aguacate cultivar Hass en cuatro municipios del Estado de México, durante el periodo de enero a junio de 2017. La distribución de la antracnosis se ajustó a modelos gaussianos y exponenciales en la mayoría de los casos. Los mapas de infestación realizados mediante krigeado muestran más de un centro de agregación de la enfermedad. Este análisis permitió estimar la superficie infestada: se encontró una infestación de más del 50% en los primeros muestreos y de hasta un 98% en los muestreos de junio en todas las zonas estudiadas. © 2019 Publicado por Elsevier Espana, S.L.U. en nombre de Asociacion Argentina de Microbiologıa. Este es un art´ıculo Open Access bajo la licencia CC BY-NC-ND (http://creativecommons. org/licenses/by-nc-nd/4.0/).

Abstract Persea americana is a species of great nutritional and economic importance for Mexico, however, like any other agricultural crop, it is affected by pests and diseases that limit its worldwide commercialization. The phytopathogenic fungus Colletotrichum gloeosporioides is the causative agent of anthracnose in avocado and manifests itself in the early stages of fruit development as well as in post-harvest and storage, under conditions of high relative humidity (80%) and at temperatures from 20°C, causing losses economic up to 20% of production. Applying geostatistical methods the present study aims to define the spatial distribution of anthracnose in Hass avocado fruits in four municipalities of the State of Mexico during the period from January to June 2017. The results show that the distribution of anthracnose was adjusted to gaussian and exponential models in most, the infestation maps made through the kriging show more than one centerof aggregation of the disease, based on it the infested surface was estimated, finding an infestation of more than 50% in the first samples and up to 98% in the samplings belonging to the month of June in all the areas studied. © 2019 Published by Elsevier Espana, S.L.U. on behalf of Asociación Argentina de Microbiología. This is an open access article under the CC BY-NC-ND license (https://creativecommons.org/ licenses/by-nc-nd/4.0/).

Spatial variability of precompression stress and volumetric water content of a red-yellow latosol (Oxisol) / Variabilidade espacial da pressão de preconsolidação e da umidade volumétrica de um latossolo vermelho amarelo

The intensive agricultural machinery traffic to which soils under coffee crops are exposed may cause significant changes in soil structure. The goals of this study were to: a) characterize the spatial variability of precompression stress (σp) and volumetric water content (θ) of a Red-Yellow Latosol (Oxisol) and determine their spatially dependent structures using ordinary kriging semivariograms; b) using this date to construct a load-bearing capacity (LBC) map for the site, identifying the depth with higher ; and c) use the map as a decision support tool regarding agricultural machinery logistic management in coffee plantations. The research was conducted on an experimental farm of the Agricultural Research Company of Minas Gerais (Epamig), Patrocinio, MG, on a clayey Oxisol. Samples were collected from 2.0 × 1.5 m sampling pits at intersections of a rectangular 40 × 150 m grid, totaling 28 sampling points. The soil samples were collected with metal rings (0.0254 m in height and 0.0630 m in diameter) using an Uhland sampler. Sampling at each pit was at 3 layers: 0.00­0.03 m, 0.10­0.13 m, and 0.25­0.28 m, and seven samples arranged in a matrix form were collected at each layer, totaling 588 samples. The σp and θ showed a spatially dependent structure. The depth of 0.00­0.03 m showed higher LBC, indicating that this depth was more compacted compared to other depths. Based on the σp maps, the tractor and the combine should not traffic into the areas at water content of 0.45 m3 m-3, because the soil has a LBC of 200 kPa. If this condition is not respected, additional compaction will occur. Tractor traffic is permissible when the water content reaches values less than 0.36 m3 m-3, whereas for combine traffic, the water content values must be less than 0.30 m3 m-3.

O tráfego intenso de máquinas agrícolas sobre solos cultivados com cafeeiro pode causar alterações significativas na estrutura do solo. Os objetivos deste estudo foram: a) avaliar e identificar a variabilidade espacial da pressão de preconsolidação (σp) e a umidade volumétrica (θ) de um Latossolo Vermelho-Amarelo (LVA) e determinar a estrutura de dependência espacial através de semivariogramas utilizando a krigagem ordinária; b) verificar no mapa, a profundidade com maior capacidade de suporte de carga (CSC) do solo ; c) dar suporte a tomada de decisão no manejo do maquinário agrícola da cultura cafeeira e utilizar esses mapas para a logística de manejo de trafego de máquinas. O experimento foi conduzido num LVA textura muito argilosa na Fazenda da Epamig (Empresa de Pesquisa Agropecuária de Minas Gerais), em Patrocínio-MG. A amostragem foi realizada numa malha retangular (150 m x 40 m), onde foram coletadas amostras indeformadas em anéis metálicos de 0,0254 m de altura por 0,0630 m de diâmetro com auxílio do amostrador tipo Uhland em 28 trincheiras. Cada trincheira possui três degraus que correspondem às profundidades de (0.00-0.03; 0.10-0.13 e 0.25-0.28 m) e as dimensões dos degraus foram de 2,0 x 1,5 m, coletaram-se sete amostras (disposta em forma matricial) para cada degrau, totalizando 588 amostras. A σp e θ apresentaram estrutura de dependência espacial. A profundidade de 0.00-0.03m apresentou maior CSC indicando que esta profundidade está mais compactada em relação as demais profundidades. Com base nos mapas da σp, o trator e a colhedora não devem trafegar na área para umidade volumétrica igual a 0.45 m3 m-3, pois o solo possui CSC de 200 kPa. Se essa condição for desrespeitada, a compactação adicional poderá ocorrer. Para o tráfego do trator, deve-se esperar por um período, para que a θ atinja um valor menor que 0.36 m3 m-3, enquanto que para o tráfego da colhedora deve-se esperar a θ atingir valor menor que 0.30 m3 m-3.

Soil chemical attributes variability under tillage and no-tillage in a long-term experiment in southern Brazi l / Variabilidade dos atributos químicos do solo em sistemas de plantio distintos em experimento de longa duração no sul do Brasil

Spatial variability of chemical attributes in soils tends to increase under the no-tillage system. However, data regarding this variability under this system for a long time are scarce. This study aimed to evaluate the spatial variability of phosphorus (P), organic matter (OM), pH, aluminum (Al), calcium (Ca), magnesium (Mg), potassium (K), base saturation (V%), sulfur (SO42-), and micronutrients: boron (B), iron (Fe), copper (Cu), manganese (Mn), and zinc (Zn) in soil with conventional tillage and no-tillage treatments for more than three decades in order to optimize fertilization and maintain soil conservation. Samples of Rhodic Hapludox were collected at 0.00­0.20 m depth. Sampling points were spaced 5 m apart, totaling 100 points. All sampling points were georeferenced according to the tillage system. Samples were analyzed chemically. Semivariograms and kriging were used in descriptive statistics and geostatistics analyses of the results. Only micronutrients and P presented distinct spatial variability between the treatments. Most of the analyzed chemical attributes presented spatial dependency in both tillage systems, which may help determine the best form of sampling and indicate adequate fertilization. According to semivariogram range values of OM, Mg, and K, the same sampling strategy may be adopted for analysis in both tillage systems, yet for Al, Ca, and pH the strategy should be specific

A variabilidade espacial dos atributos químicos tende a aumentar em solos sob sistema plantio direto. No entanto faltam informações sobre essa variabilidade em solos submetidos a esse sistema por longo período.Este estudo objetivou avaliar a variabilidade espacial do fósforo (P), matéria orgânica (MO), pH, alumínio (Al), cálcio (Ca), magnésio (Mg), potássio (K), saturação por bases (V%), enxofre (SO42-) e dos micronutrientes boro (B), ferro (Fe), cobre (Cu), manganês (Mn) e zinco (Zn) em plantio direto e plantio convencional adotados há mais de três décadas, como forma de otimizar a aplicação de fertilizantes e preservar a adoção do manejo conservacionista. Amostras de Latossolo Vermelho distrófico coletadas na profundidade de 0,00 - 0,20 m, e 5 m de distância entre os 100 pontos amostrados e georreferenciados por manejo, foramanalisadas quimicamente e os resultados submetidos à estatística descritiva e à geoestatística, através do estudo de semivariogramas e krigagem. Apenas os micronutrientes e o P apresentaram variabilidade espacial distinta entre os manejos. A maioria dos atributos químicos estudados apresentou dependência espacial em ambos os sistemas de manejo, o que pode determinar a melhor forma de amostragem e indicar a adubação adequada. De acordo com os alcances atribuídos aos semivariogramas, para MO, Mg e K pode-se adotar a mesma estratégia de amostragem para análise em ambos sistemas de manejo e para Al, Ca e pH deve ser específica

Spatial variation of chemical attributes in archaeological dark earth under cocoa cultivation in Western Amazon / Variação espacial de atributos químicos em terra pr eta arqueológica sob cultivo de cacau na amazonia ocide ntal

As Terras Pretas Arqueológicas (TPAs) são solos de fertilidade que notoriamente são superiores à grande maioria dos solos típicos da região amazônica. Assim, o uso das TPAs intensificou-se devido às características químicas apresentadas por esses solos, como a alta concentração de nutrientes (fósforo, cálcio, magnésio). Nesse sentido, o objetivo deste estudo foi avaliar a distribuição espacial dos atributos químicos do solo em uma área de Terra Preta Arqueológica sob cultivo de cacau no município de Apuí (AM). O mapeamento de uma malha de 42 x 88 m, com espaçamento irregular de 6 x 8 m, totalizando 88 pontos, foi realizado e, em seguida, foram coletadas amostras de solo nas profundidades de 0,0-0,05; 0,05-0,10; (pH, O.C, Sto. C, (H + Al), P, K, Ca, Mg, SB, CEC e V%). Os dados foram analisados por meio de técnicas estatísticas descritivas e geoestatísticas. Os valores médios e medianos foram ajustados aos valores próximos, indicando distribuição normal, enquanto os atributos químicos do solo foram ajustados aos modelos de semivariogramas esféricos e exponenciais. A maioria dos atributos apresentou coeficiente de variação (CV) entre 12,1 e 60%, caracterizada como variabilidade média, as variáveis do estudo apresentaram diferentes faixas e a maioria teve forte dependência espacial. As técnicas geoestatísticas utilizadas permitiram os ajustes dos modelos teóricos que melhor representaram a semivariância experimental, possibilitando a construção de mapas temáticos da distribuição espacial dos valores dos atributos da área estudada.

The Black Archaeological Lands (TPAs) are fertility soils that are notoriously superior to the vast majority of typical soils in the Amazon region. Thus, the use of TPAs ​​has intensified due to the chemical characteristics presented by these soils, such as the high concentration of nutrients (phosphorus, calcium, magnesium). In this sense, the objective of this study was to evaluate the spatial distribution of the chemical attributes of the soil in an area of ​​Terra Preta Arqueológica under cocoa cultivation in the city of Apuí (AM). The mapping of a mesh of 42 x 88 m, with irregular spacing of 6 x 8 m, totaling 88 points, was carried out and, then, soil samples were collected at depths of 0.0-0.05; 0.05-0.10; (pH, O.C, Sto. C, (H + Al), P, K, Ca, Mg, SB, CEC and V%). The data were analyzed using descriptive and geostatistical statistical techniques. The mean and median values ​​were adjusted to the close values, indicating normal distribution, while the chemical attributes of the soil were adjusted to the models of spherical and exponential semivariograms. Most of the attributes presented a variation coefficient (CV) between 12.1 and 60%, characterized as average variability, the study variables had different ranges and most had a strong spatial dependence. The geostatistical techniques used allowed the adjustments of the theoretical models that best represented the experimental semivariance, allowing the construction of thematic maps of the spatial distribution of the values ​​of the attributes of the studied area.

Variability and spatial correlation of aggregates and organic carbon in indian dark earth in Apuí region, AM / Variabilidade e correlação espacial de agregados e carbono orgânico em terra preta de índio na região do Apuí, AM

Having in mind the importance of knowing the variability and spatial correlation of soil properties in Indian Dark Earth (IDE), we evaluated in this study the variability and the spatial correlation of aggregates and carbon in an Ultisol under coffee cultivation in southern Amazonas. It was established a 48 x 88 m sampling grid spaced 06 x 08 m, totalling 88 sampling points. Then soil samples were collected at: 0.0-.05, 0.05-0.10, and 0.10-0.20 m layers. The spatial variability of the Mean Weighted Diameter (MWD) attributes, aggregates > 2 mm, < 2 mm, bulk density (BD) and organic carbon (OC) was analyzed by adjusting the simple semivariograms, while spatial correlations of the OC with aggregates and BD were analyzed by cross-semivariogram. We could conclude that there was spatial dependence in the variables, wherein the largest ones were observed at 0.0-0.05 m, except for Mean Weighted Diameter (MWD) and aggregates greater than 2.00 mm with larger range of values in depth from 0.05-0.10 and 0.10-0.20 m. The mean weight diameter and aggregate class attributes greater than 2.00 mm had negative spatial correlation with organic carbon at 0.0-0.05 m, while the smaller aggregates than 2.00 mm classes and bulk density correlated positively with organic carbon at 0.0-0.05 m and 0.10-0.20 m.

Considerando a importância do conhecimento da variabilidade e correlação espacial dos atributos do solo em Terra Preta de Índio (TPIs), avaliou-se neste trabalho a variabilidade e a correlação espacial de agregados e carbono em um Argissolo Amarelo eutrófico sob cultivo de café na região sul do Amazonas. Foi estabelecido um grid amostral com dimensões de 48 x 88 m e espaçamentos de 06 x 08 m, totalizando 88 pontos amostrais. Em seguida, foram coletadas amostras de solos nas profundidades: 0,0-0,05, 0,05-0,10, e 0,10-0,20 m. A variabilidade espacial dos atributos diâmetro médio ponderado (DMP), agregados > 2 mm, < 2 mm, densidade do solo (Ds) e carbono orgânico (CO) foi analisada ajustando os semivariogramas simples, enquanto as correlações espaciais do CO com agregados e Ds foram analisadas por semivariogramas cruzados. Concluiu-se que as variáveis apresentaram dependência espacial, e os maiores alcances são constatados na profundidade 0,0-0,05 m, exceto para DMP e agregados maiores que 2,00 mm com maiores valores de alcance na profundidade 0,05-0,10 e 0,10-0,20 m. Os atributos diâmetro médio ponderado e classes de agregados maior que 2,00 mm apresentam correlação espacial negativa com carbono orgânico na profundidade 0,0-0,05 m, enquanto a classes de agregados menor que 2,00 mm e densidade do solo apresentam correlação positiva com carbono orgânico nas profundidades 0,0-0,05 m e 0,10-0,20 m.

Space variability of phenological indicators of common bean crop / Variabilidade espacial dos índices fenológicos da cultura do feijão

Common bean (Phaseolus vulgaris L.) has a representative agricultural holding, not only for the economic value of its production, but also for the large area of growing in Brazil. In the harvest 2016/17, this work was conducted in a Quartzarenic Neosol in the municipality of Cassilândia, MS. The objective of this work was to characterize the structure and magnitude of the spatial distribution of phenological indices of the common bean crop and to map the phenological indices in order to visualize the spatial distribution and to evaluate the spatial correlation among common bean yield and plant variables: grain yield (YIE), mass of one hundred grains (MHG), number of grains per plant (NG), number of grains per pod (NGP), number of pods per plant (NP), dry matter (DM), plant length (PL) and stem diameter (SD), sampled in a grid of 117 georeferenced points (81 points of base grid and 36 points of higher density grid). Analysis of these data through statistical and geostatistical techniques made it possible to verify that the production and yield components presented spatial dependence. There was a positive spatial correlation among common bean yield and the mass of one hundred grains, number of grains per pod and plant length, demonstrating that they have a strong spatial dependence.

A cultura do feijoeiro (Phaseolus vulgaris L.) tem representativa exploração agrícola, não só pelo valor econômico de sua produção, como também pela grande área de cultivo no Brasil. No ano agrícola de 2016/17, este trabalho foi conduzido em um Neossolo Quartzarênico no município de Cassilândia, MS. O trabalho objetivou caracterizar a estrutura e a magnitude da distribuição espacial de índices fenológicos da planta em lavoura de feijão e realizar o mapeamento desses índices fenológicos, de forma a visualizar a distribuição espacial, e avaliar a correlação espacial existente entre a produtividade do feijoeiro e as variáveis da planta: produtividade de grãos (PG), massa de cem grãos (MC), número de grãos por uma planta (GP), número de grãos por uma vagem (GV), número de vagem por uma planta (VP), massa seca de uma planta (MS), comprimento da planta (CO), diâmetro do colmo (DC), amostrados em uma malha de 117 pontos georreferenciados (81 pontos da malha base e 36 pontos de malha com maior densidade). A análise destes dados por meio das técnicas estatísticas e da geoestatística possibilitaram constatar que os componentes de produção e produtividade do feijão apresentaram dependência espacial. Houve destaque na correlação espacial positiva entre a produtividade do feijoeiro e a massa de cem grãos, grãos por vagem e comprimento da planta, demonstrando que as mesmas possuem uma dependência espacial forte.

Some aspects about the spatial dependence index for variability of soil attributes / Alguns aspectos sobre o índice de dependência espacial para variabilidade de atributos do solo

ABSTRACT: The main purpose of this article was to evaluate the behavior and relationship of the range and components of SDI (Spatial Dependence Index) in general and in function of field factors such as soil types, type of attribute and soil layers. This evaluation was based on real data collected in national journals. It was noticed that the parameter range, in general and for different field factors, presented asymmetric positive behavior. The components of the SDI showed approximately symmetrical behavior. The SDI can capture the range behavior more intensely (the spatial variability behavior in the horizontal direction of the semivariogram), and, in a less intense way, the behavior of the contribution and sill parameters (the spatial dependence behavior in the vertical direction of the semivariogram). Thus, the SDI describes the behavior of spatial dependence of the total set of aspects of the semivariogram.

RESUMO: O objetivo deste artigo foi avaliar o comportamento e o relacionamento do alcance e dos componentes do IDE (Índice de Dependência Espacial) em geral e em função de fatores de campo, tais como tipos de solo, tipo de atributo e profundidades de solo. Esta avaliação foi baseada em dados reais coletados em periódicos nacionais. Foi observado que o parâmetro alcance, em geral, e para diferentes fatores de campo, apresentou comportamento de assimetria positiva. Os componentes do IDE apresentaram comportamento aproximadamente simétrico. O IDE consegue captar, de forma mais intensa, o comportamento do alcance (o comportamento da variabilidade espacial no sentido horizontal do semivariograma), e, de forma menos intensa, o comportamento dos parâmetros contribuição e patamar (o comportamento da dependência espacial no sentido vertical do semivariograma). Assim, o IDE descreve o comportamento de dependência espacial de toda a gama de aspectos do semivariograma.

Object-based change detection using semivariogram indices derived from NDVI images: The environmental disaster in Mariana, Brazil / Detecção de mudanças baseada em objetos utilizando indices do semivariograma derivados de imagens NDVI: O disastre ambiental em Mariana, Brasil

ABSTRACT Object-based change detection is a powerful analysis tool for remote sensing data, but few studies consider the potential of temporal semivariogram indices for mapping land-cover changes using object-based approaches. In this study, we explored and evaluated the performance of semivariogram indices calculated from remote sensing imagery, using the Normalized Differential Vegetation Index (NDVI) to detect changes in spatial features related to land cover caused by a disastrous 2015 dam failure in Brazil's Mariana district. We calculated the NDVI from Landsat 8 images acquired before and after the disaster, then created objects by multiresolution segmentation analysis based on post-disaster images. Experimental semivariograms were computed within the image objects and semivariogram indices were calculated and selected by principal component analysis. We used the selected indices as input data to a support vector machine algorithm for classifying change and no-change classes. The selected semivariogram indices showed their effectiveness as input data for object-based change detection analysis, producing highly accurate maps of areas affected by post-dam-failure flooding in the region. This approach can be used in many other contexts for rapid and accurate assessment of such land-cover changes.

RESUMO Recentemente, variáveis geoestatísticas derivadas de imagens de sensoriamento remoto ganharam espaço dentre os procedimentos de detecção de mudanças, porém, o potencial temporal destas variáveis para o mapeamento das mudanças baseado na análise por objetos ainda é pouco estudado. Neste estudo, o desempenho de um conjunto de índices calculados de semivariogramas derivados de imagens NDVI bitemporais para detectar mudanças na cobertura do solo foi analisado e avaliado. O município de Mariana foi selecionado para teste e validação da metodologia devido ao grande impacto ocasionado pelo desastre. O processo iniciou-se com a aquisição de imagens Landsat 8 antes e após o desastre e o cálculo do NDVI. Os objetos foram criados através da segmentação em multiresolução baseada na imagem pós-desastre. Os semivariogramas experimentais foram gerados dentro de cada objeto e os índices foram extraídos e selecionados através da análise de componentes principais. Os índices selecionados foram utilizados como dados de entrada para o algoritmo support vector machines para a classificação de áreas de mudança e não mudança. Os índices selecionados se mostraram efetivos para a detecção de mudanças, indicando a possibilidade de utilização para a detecção de mudanças baseada em objetos, resultando em um mapa precisos das áreas inundadas afetadas pelo desastre. Esta abordagem pode ser usada em muitos outros contextos para uma avaliação rápida e precisa de tais mudanças na cobertura do solo.

Estatística descritiva e estacionaridade em variáveis geoquímicas ambientais / Descriptive statistics and stationarity in environmental geochemical variables

RESUMO O objetivo deste trabalho foi analisar os dados de atributos geoquímicos a fim de verificar sua estacionaridade e correlacionar a normalidade estatística com o uso da técnica de krigagem ordinária. A escolha da krigagem ordinária como método geoestatístico aplicado ao trabalho deve-se ao fato de essa ser aconselhada para a realização de estudos em áreas onde existam dados com variáveis que possam apresentar dependência espacial, como é o caso das variáveis geoquímicas, e por ser indicada para dados que apresentam estacionaridade. A metodologia utilizada para a realização desta pesquisa envolveu, além da revisão de literatura, a obtenção de dados dos metais-traço (Cu, Zn, Mn, Fe, Cr e Pb) extraídos parcialmente de amostras superficiais (0 a 10 cm) de solos e sedimentos coletados em campo. Também foram determinados os valores de pH, salinidade, nitrogênio total, fósforo, matéria orgânica e granulometria. Foram conduzidas análises estatísticas, construções de semivariogramas, aplicação da krigagem ordinária e, por fim, validação cruzada para medir a incerteza da medição prévia dos dados. Neste trabalho, por meio dos variogramas, comprovou-se que, apesar de os dados não serem normais, eles apresentaram estacionaridade. Além disso, o parâmetro da estatística descritiva que mais possui correlação direta com a krigagem ordinária é a variância.

ABSTRACT The aim of this work was to analyze geochemical data in order to check their stationarity and to correlate the statistical normality using the ordinary kriging technique. The ordinary kriging technique was chosen as the geostatistical method applied to work because such technique is advised for studies in areas where there are data with variables that might present spatial dependence, like the geochemical variables, and also because it is indicated for data presenting stationarity. The methodology used for this research involved, besides literature review, data collection of trace metals (Cu, Zn, Mn, Fe, Cr and Pb) that were partially extracted from surface samples (0 to 10 cm) of soils and sediments collected in the field. We also determined the values of pH, salinity, total nitrogen, phosphorus, organic matter and particle size. Statistical analyzes, semivariogram development, ordinary kriging use and, lastly, cross validation were performed to measure the uncertainty of the previous measurement of data. It was found, in this work, by means of the variograms that although data were ordinary, they showed stationarity. In addition, the parameter of descriptive statistics that mostly correlates directly with the ordinary kriging is variance.

Hybrid kriging methods for interpolating sparse river bathymetry point data / Métodos híbridos de krigagem para interpolação de levantamentos batimétricos fluviais

ABSTRACT Terrain models that represent riverbed topography are used for analyzing geomorphologic changes, calculating water storage capacity, and making hydrologic simulations. These models are generated by interpolating bathymetry points. River bathymetry is usually surveyed through cross-sections, which may lead to a sparse sampling pattern. Hybrid kriging methods, such as regression kriging (RK) and co-kriging (CK) employ the correlation with auxiliary predictors, as well as inter-variable correlation, to improve the predictions of the target variable. In this study, we use the orthogonal distance of a (x, y) point to the river centerline as a covariate for RK and CK. Given that riverbed elevation variability is abrupt transversely to the flow direction, it is expected that the greater the Euclidean distance of a point to the thalweg, the greater the bed elevation will be. The aim of this study was to evaluate if the use of the proposed covariate improves the spatial prediction of riverbed topography. In order to asses such premise, we perform an external validation. Transversal cross-sections are used to make the spatial predictions, and the point data surveyed between sections are used for testing. We compare the results from CK and RK to the ones obtained from ordinary kriging (OK). The validation indicates that RK yields the lowest RMSE among the interpolators. RK predictions represent the thalweg between cross-sections, whereas the other methods under-predict the river thalweg depth. Therefore, we conclude that RK provides a simple approach for enhancing the quality of the spatial prediction from sparse bathymetry data.

RESUMO Modelos de terreno de rios são usados para análise de mudanças geomorfológicas e para simulações hidrológicas. Estes modelos são interpolados a partir de pontos batimétricos. A batimetria fluvial é geralmente conduzida através de seções transversais, o que pode acarretar em uma malha amostral esparsa. Métodos híbridos de krigagem, como krigagem por regressão (KR) e co-krigagem (CK), empregam a correlação com preditores auxiliares, além da auto-correlação entre variáveis, na predição da variável resposta. Neste estudo, sugere-se que a distância ortogonal de um ponto até a linha de centro do talvegue de um rio pode ser usada como covariável para KR e CK. Considerando-se que a variabilidade da cota do leito do rio é abrupta transversalmente a direção do fluxo, espera-se que quanto maior a distância euclidiana de um ponto até o talvegue, maior será sua elevação. O objetivo deste estudo foi avaliar o uso da covariável proposta em métodos híbridos de krigagem para a predição espacial da topografia do leito de rios. Para tanto, foi realizada uma validação externa, em que seções transversais foram usadas para interpolação e dados levantados entre as seções consistiram na amostra de teste. Os resultados da KR e CK foram comparados aos da krigagem ordinária. A KR apresentou a menor REQM. No mapa resultante da KR, o talvegue foi preservado nas lacunas não amostradas entre as seções, enquanto os demais métodos subestimaram a profundidade do talvegue nestes espaços. Assim, conclui-se que a KR pode melhorar a predição espacial de dados batimétricos fluviais.

Spatial and temporal distribution of species diversity in semideciduous seasonal forests with occurrence of fire / Distribución espacial y temporal de la diversidad de especies en un bosque estacional semideciduo con un episodio de incendio

AbstractThe spatial and temporal distributions of vegetation can be influenced by edaphic and environmental factors as well as disturbances. The present study was undertaken to characterize the distribution and spatial dependence of species diversity in a fragment of semideciduous seasonal forest after a disturbance event by fire and to identify changes during natural regeneration. The fire, occurred in 1999, and surveys were undertaken in 2000, 2005 and 2014 in 26 geo-referenced plots (400 m²) distributed along two transects. The Fire Damage Index (FDI) in each plot was based on a scale from 0 to 5, and all of the arboreal individuals with diameters at breast height (DBH) ≥ 5 cm were identified in each plot. Species diversity was calculated using the Shannon index (H'), and species richness (S) was estimated using species accumulation curves; variations between the surveys compared based on the means and standard errors. The S, H' and FDI data were submitted to non-spatial exploratory and geostatistical analyses. After adjusting the semivariograms, the levels of spatial dependence were classified, and interpolation of the variables values were performed using ordinary kriging to characterize their spatial distributions in the form of maps. Spatial analysis was used to identify and characterize differences between the post-fire surveys in terms of the configurations of the arboreal community. The fire event influenced the spatial and temporal structures of the variables S and H'. These variables showed spatial dependence and aggregated distributions, with reduction in the distance under spatial influences and a uniformity of individuals distribution in the forest fragment at the different surveys. This research characterized the distribution and spatial dependence of the variables S and H' in a forest fragment after a fire event and the alterations in the arboreal community structure during natural regeneration.

ResumenLa distribución espacial y temporal de la vegetación esta influenciada por factores edáficos, ambientales y alteraciones. En ese contexto, el objetivo del estudio fue caracterizar el patrón de distribución y la dependencia espacial de la diversidad de especies en un fragmento de bosque estacional semideciduo, donde ocurrió un incendio, e identificar alteraciones a lo largo del proceso de regeneración natural. El incendio ocurrió en 1999, y en 2000, 2005 y 2014, se hicieron inventarios florísticos en 26 parcelas (400 m2) georreferenciadas y distribuidas en dos transectos. El índice de daño del fuego (IDF) se evaluó en cada parcela con base en una escala de 0 a 5. En cada inventario se identificaron las plantas arbóreas con diámetro a altura de pecho (DAP) ≥ 5 cm, calculando la diversidad de especies con base en el índice de Shannon (H') y estimando la riqueza de especies (S). Esas variaciones fueron analizadas utilizando la curva de acumulación de especies, y la variación entre inventarios fue comparada con base en las medias y el error estándar. Los datos de S, H' e IDF fueron sometidos a análisis exploratorio no espacial y a análisis geoestadístico. Después del ajuste de los semivariogramas, se clasificó el nivel de dependencia espacial y se hizo la interpolación de datos por Krigeado ordinario para caracterizar, en mapas, el patrón de distribución espacial de esas variables. Se encontró diferencia entre los inventarios realizados después del incendio, como en la configuración de la comunidad arbórea, por lo que el análisis espacial permitió identificar y caracterizar este hecho. El incendio afectó la estructura espacial y temporal de las variables S y H'. Estas variables presentan dependencia espacial y la distribución agregada, con distancias más pequeñas bajo influencia espacial, y la uniformidad de la distribución de los individuos en el fragmento de bosque durante las distintas etapas. Por último, se caracterizó la distribución y la dependencia espacial de las variables S e H 'en un fragmento de bosque después de un incendio, y las alteraciones en la estructura de la comunidad arbórea durante la regeneración natural.

Spatial assessment of potential ecological risk of heavy metals in soils from informal e-waste recycling in Ghana

The rapidly increasing annual global volume of e-waste, and of its inherently valuable fraction, has created an opportunity for individuals in Agbogbloshie, Accra, Ghana to make a living by using unconventional, uncontrolled, primitive and crude procedures to recycle and recover valuable metals from this waste. The current form of recycling procedures releases hazardous fractions, such as heavy metals, into the soil, posing a significant risk to the environment and human health. Using a handheld global positioning system, 132 soil samples based on 100 m grid intervals were collected and analysed for cadmium (Cd), chromium (Cr), copper (Cu), mercury (Hg), lead (Pb) and zinc (Zn). Using geostatistical techniques and sediment quality guidelines, this research seeks to assess the potential risk these heavy metals posed to the proposed Korle Ecological Restoration Zone by informal e-waste processing site in Agbogbloshie, Accra, Ghana. Analysis of heavy metals revealed concentrations exceeded the regulatory limits of both Dutch and Canadian soil quality and guidance values, and that the ecological risk posed by the heavy metals extended beyond the main burning and dismantling sites of the informal recyclers to the school, residential, recreational, clinic, farm and worship areas. The heavy metals Cr, Cu, Pb and Zn had normal distribution, spatial variability, and spatial autocorrelation. Further analysis revealed the decreasing order of toxicity, Hg>Cd>Pb> Cu>Zn>Cr, of contributing significantly to the potential ecological risk in the study area.

DISTRIBUCIÓN ESPACIAL, COMPOSICIÓN Y DENSIDAD DE EDAFOFAUNA EN HOJARASCA DE BOSQUE Y CAFETAL (MONTENEGRO, COLOMBIA) / Spatial Distribution, Composition and Density of the Leaflitter Edaphofauna in Forest and Coffee Plantation (Montenegro, Colombia)

Se determinaron los cambios en composición, densidad y distribución espacial de la edafofauna de hojarasca entre dos usos de suelo (cafetal y bosque subandino) para entender el efecto del uso de suelo. Asimismo, se establecieron las relaciones entre las variables de la edafofauna y atributos de la hojarasca. Se colectaron muestras de hojarasca superficial durante la época seca y se extrajo la edafofauna, manualmente y mediante embudos Berlese modificados. Esta se identificó hasta el menor nivel taxonómico posible y se clasificó en grupos funcionales. La composición faunística fue diferente entre los usos de suelo. La similitud en la composición fue menor al 50 % en dos de los tres grupos funcionales encontrados, pero la densidad total no presentó diferencias. En general hubo coeficientes de variación altos en la densidad faunística intermuestral para los diferentes grupos taxonómicos y, en consecuencia, una baja detectabilidad del patrón espacial. Sin embargo, se detectó un patrón espacial en la abundancia para Coleoptera, larvas de Diptera y biomasa de hojarasca, en el cafetal. La estructura fue más difícil de detectar en bosque, por los pequeños parches en este. La varianza de la abundancia fue mayor en el cafetal, para todos los grupos, excepto Acaridida y Formicidae. La densidad de Formicidae aumentó al hacerlo la humedad de la hojarasca, pero solo en cafetal. Se concluye que el establecimiento del cafetal afecta de manera diferencial a taxones y grupos funcionales de la edafofauna. Esto sienta las bases para definir los efectos del cambio de uso de suelo y los factores ambientales asociados con la distribución espacial de la edafofauna.

Changes in composition, density and spatial distribution on the litter edaphic fauna between two land uses (a coffee plantation and a tropical Andean forest) were established, in order to understand the soil use effects. Also, the relationships between the biotic variables and some litter features were evaluated. Superficial litter samples were collected during dry season. The edaphic fauna was extracted manually and with modified Berlese funnels; it was identified up to the lower taxonomic level and classified in functional groups. Differences in fauna composition were observed between the two land uses. Composition similarity was lower than 50 % for almost all functional groups found, but total density did not differ between the land uses. In general, there was a high variation in density among samples resulting in higher coefficients of variation for different taxonomic groups and consequently low detectability of spatial pattern. However a spatial pattern was detected only for Coleoptera, Diptera larvae and litter biomass in the coffee plantation. The structure in forest was more difficult to detect, due to the smaller patches on it. The variance of abundance was significantly higher in the coffee plantation for all groups, except Acaridida and Formicidae. Density of Formicidae increased with the enhancement of litter moisture, but only in the coffee plantation. In conclusion, the establishment of coffee plantation differentially affects taxa and functional groups of soil fauna. This is the basis for establishing the effects of soil change use, and the environmental factors associated with spatial distribution of edaphic fauna.